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Abstract. A 4-dimensional relativistic positioning system for a general 
spacctime is constructed by using the so called "emission coordinates." The results 
apply in a small region around the world line of an accelerated observer carrying 
a Fermi triad, as described by the Fermi metric. In the case of a Schwarzschild 
spacetime modeling the gravitational field around the Earth and an observer at 
rest at a fixed spacetime point, these coordinates realize a relativistic positioning 
system alternative to the current GPS system. The latter is indeed essentially 
conceived as Newtonian, so that it necessarily needs taking into account at least 
the most important relativistic effects through Post-Newtonian corrections to 
work properly. Previous results concerning emission coordinates in flat spacetime 
are thus extended to this more general situation. Furthermore, the mapping 
between spacetime coordinates and emission coordinates is completely determined 
by means of the world function, which in the case of a Fermi metric can be 
explicitly obtained. 



PACS number: 04.20.Cv 

1. Introduction 

Currently there is a growing interest in the construction of an emission coordinate 
system for the Earth in order to improve the current positioning systems (GPS, 
GLONASS) [2 [2]. In fact, the latter are essentially conceived as Newtonian, hence 
based on a classical (Euclidean) space and absolute time, over which some relativistic 
corrections are added via the Post-Newtonian formalism. On the contrary, emission 
coordinates allow the definition of relativistic positioning systems whose study based 
on the framework and the concepts of general relativity was initiated by Coll and 
collaborators several years ago [3J HI [U GO El [5] . 

At present global positioning systems consist of a non-inertial spatial reference 
frame for navigation that co-rotates with the Earth and is geocentric (the ECEF, 
Earth Centered Earth Fixed system), and on the coordinate time of a local inertial 
"star-fixed" reference frame whose origin lies at the Earth's center of mass and which is 
freely falling with it (the ECI, Earth Centered Inertial system). Although clock speeds 
are small in comparison with the speed of light and gravitational fields are weak near 
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the Earth, they give rise to significant relativistic effects. The most important ones 
affecting the rate of clocks (to the order 1/c 2 , which is just the order of approximation 
used in GPS) are first and second order Doppler frequency shifts of clocks due to their 
relative motion, gravitational frequency shifts due to the Earth's mass, and the Sagnac 
effect due to the Earth's rotation [2 [2 [9]. If these corrections were not made, the 
system would not be operational after a few minutes. In a day of operation, they would 
produce an error of more than 11 km in the horizontal positioning of the receiver. In 
a week, the error accumulated in the vertical positioning would be approximately 5 
km. Therefore, fully-relativistic positioning systems are not only conceptually simpler 
than GPS systems, but also more accurate, because no corrections are needed at all, 
whereas Post-Newtonian corrections to current positioning systems are included up to 
the order 1/c 2 only. 

Recently the project SYPOR [3] (SYsteme de POsitionnement Relativiste) has 
proposed to endow the constellation of satellites of GALILEO with the necessary 
elements to constitute by itself a primary, autonomous positioning system for the Earth 
and its neighbors, referring to the capability of the constellation to provide complete 
relativistic metric information, i.e., to describe both the kinematics and the dynamics 
of the constellation itself as well as of the users. In this primary positioning system, 
an observer at any event in a given spacetime region can know its proper coordinates. 
The system is also autonomous or autolocated if any receiver determines its spacetime 
path as well as the trajectories of the satellites solely on the basis of the information 
received during a proper time interval. Hence, a constellation of satellites with clocks 
that interchange their proper time among themselves and with Earth receivers is a 
fully relativistic system. According to the SYPOR project the GALILEO system 
would not need "relativistic" corrections. Giving a theoretical contribution to that 
project is the main motivation of the present work. 

The simplest relativistic positioning system is the one formed by electromagnetic 
signals broadcasting the proper times t a (A = 1, ... ,4) of four independent clocks 
carried by satellites [J which move along geodesic (i.e., freely falling) world lines. The 
above signals, parameterized by the proper time of the clocks, realize four emission 
coordinates {t a }. More precisely, let us consider an arbitrary spacetime covered by 
the coordinate patch {x a }. Let P be a generic spacetime event with coordinates X a 
and let Pa be a generic point with coordinates X% along the world line of the Ath 
satellite. The condition ensuring that the signals emitted by the four satellites at Pa 
meet the receiver location at P is given by 

n(x A ,x) = o, x° A <x°, (i.i) 

where Q(X A , X) is the world function connecting the receiver and emission points and 
P is in the strict causal future of Pa- For any given background metric, the world 
function is defined as half the square of the spacetime distance between two generic 
points xa and xb connected by a geodesic path (see Eq. (1) of Chap. 2 in Ref. [10]) 



1 f 1 

n(xA,x B ) = - J 



dx 1 " d/„ , „. 

dA , (1.2) 



where x a (X) satisfies the geodesic equation and the affine parameter is such that 
x a (0) = x^ and x a (l) = x B . Eq. Ql.lj) is a system of four equations which must 
be solved for the unknown coordinates X a of P in terms of the satellite coordinates 

X Hereafter we will always use the single word "satellite" to mean "clock carried by satellite," for 
convenience. 
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X A , supposed to be known in terms of the proper times t a of the satellites, i.e., 
the emission coordinates of the point P. The calculation of the world function in a 
generic spacetime is not a trivial task. It is generally performed perturbatively, unless 
the solution of the geodesic equations is explicitly known, which is usually not the 
case. 

The problem of setting up such an emission coordinate system in the case of flat 
spacetime has been addressed by Coll and collaborators [6l[7l[8], with special interest in 
the 2-dimensional case, which is the simplest situation to deal with. The 4-dimensional 
case has been considered by Rovelli [TT] , who has outlined a procedure to construct a 
system of emission coordinates (introduced there with the name GPS coordinates) for a 
particular linear configuration of emitters in flat spacetime, consisting of four satellites 
moving away from the origin in different directions but at a common speed. The 
geometrical interpretation of Rovelli's construction has been discussed by Blagojevic 
et al. P] . 

In the present paper we explicitly construct emission coordinates for a general 
spacetime, in a small region around the world line of an accelerated observer carrying 
a Fermi triad, as described by the Fermi metric [13] . In particular, we study the case 
of the Schwarzschild spacetime modeling the gravitational field around the Earth and 
an observer at rest at a fixed point. 

The procedure is first outlined in flat spacetime, with a convenient choice of 
satellite motion, leading to simple explicit expressions for the metric components in 
terms of the new coordinates. This analysis is then repeated for the more interesting 
case of the Fermi metric. Emission coordinates as well as the components of the 
transformed metric are obtained as corrections to the flat spacetime ones. Since the 
constructing procedure of emission coordinates is completely general, the calculations 
can be easily extended to different choices of satellite motion. 

2. Flat spacetime 

Let us briefly review the standard construction of GPS coordinates in flat spacetime 
[11] , whose generalization to the case of Fermi background metric will be discussed in 
the next section^. 

Consider Minkowski spacetime in standard Cartesian coordinates (t, x, y, z). The 
corresponding line element is given by 

ds 2 = r! a(j dx a dx () = -dt 2 + dx 2 + dy 2 + dz 2 . (2.1) 

Let the four satellites be represented by test particles in geodesic motion. With this 
choice of coordinates timelike geodesies are straight lines 

x a A {T A )^S a A = U%T A + S% A , A=l,...,4 (2.2) 

where 

Ua = lA[dt + VAn A di] = cosh o^c^ + smhctAn l A di (2.3) 

are their (constant) 4- velocities and t a is the proper time parametrization along each 
world line. In Eq. (|2.3p ja is the Lorentz factor and the linear velocities va are related 
to the rapidity parameters a a by va = tanha^; ua denote the spacelike unit vectors 
along the spatial directions of motion. Without any loss of generality, we assume that 

§ Note that the signature conventions adopted here are different from those of 1111 . 



Emission vs Fermi coordinates 4 

the satellites all start moving from the origin of the coordinate system O; so hereafter 
we set Sq A = 0, and hence 

SI = U a A r A . (2.4) 

Let us consider now a generic spacetime point P with coordinates W a and 
the generic point Pa with coordinates S% along the world line of the Ath satellite 
corresponding to an elapsed amount of proper time t a . A photon emitted at Pa 
follows a null geodesic, i.e., the straight line 

x a (X) = W a — K a \ + S A , (2.5) 

where A is an affine parameter. Such a photon will reach P at a certain value A 
according to 

W a = K a \ + S A : , (2.6) 

implying that 

U%t a -W a = -K a X . (2.7) 

Taking the norm of both sides we get 

- {t a ) 2 + \\W\\ 2 - 2t a (U a ■ W) = , (2.8) 

since K is a null vector. 

• Emission vs spacetime coordinates 

Solving for t a and selecting the solution corresponding to the past light cone leads to 
(see Eq. (15) of Ref. [11]) 

t a = -(U A ■ W) - ^(U A ■ Wf + \\W\\ 2 . (2.9) 

These equations give the four proper times t a associated with each satellite in 
terms of the Cartesian coordinates of the generic point P in the spacetime, i.e., 
t a = T A (W a , . . . , W 3 ). The construction of emission coordinates is briefly sketched 
in Fig. 1. 

Using Eq. (|2.9[) . one can evaluate the inverse of the transformed metric 

9 AB = - ^(dr^(dr% = Ar A ■ &r* , (2.10) 

where the dual frame (d.T A ) a — dr A /dW a also satisfies the following properties 

(dr A ) a W a = t a , (dr A ) a U2 = 1 . (2.11) 
Similarly one can introduce the frame vectors 

"V *»" <.**).(»)•-#. (2.2, 



dr A J dr A ' \ dr l 

It is then easy to show [TT] [T^] that the condition g = dr A ■ dr A = is fulfilled. In 

3f ] 

W a - T A U Aa 

^ + (U A ■ W) ' 



fact, by differentiating both sides of Eq. (|2.8p with respect to W a one obtains 

(dr A ) a = ;v ! ,;. 1 . 12.1:5) 

which implies 



g AA = {dT A ndT A )a = -(r A Y + W\?-2r_ A (U A -W) = 

y y ' K ' [t a + (U A ■ W)} 2 v 1 
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The metric coefficients g AB = (d/dr A ) ■ {d/dr B ) = ^(d/dr^^d/dr^p can 
be easily obtained as well by expressing the Cartesian coordinates of P in terms of 
the emission coordinates t a , i.e., W a = W 01 ^ 1 , . . . , r 4 ). 
• Spacetime vs emission coordinates 

To accomplish this, it is enough to invert the transformation (|2.9p . However, in order 
to outline a general procedure, we start by considering the equation of the past light 
cone of the generic spacetime point P with coordinates W a given in terms of the world 
function, which in the case of flat spacetime is simply given by 

^R & t(xA,X B ) = 7jVafj{xA ~ x b)( x A ~ x b) ' ( 2 ' 15 ) 

The condition ensuring that the past light cone of P cuts the emitter world lines 
is given by 

n Rat (S A ,W) = , S° A <W°, (2.16) 

for each satellite labeled by the index A. This gives rise to a system of four quadratic 
equations in the four unknown coordinates W a of the event P of the form (|2.8|) for 
each A = 1, ... ,4. To solve this system start for example by subtracting the last 
equation from the first three equations to obtain the following system 

nfiatOSi, W) - fi flat (5 4 , W) = = -2W ■ (5, - S 4 ) - (r 1 ) 2 + (r 4 ) 2 , i = 1, 2, 3 
n fla t(S 4 , W) = = ||T?|| 2 - 2W ■ S 4 - (r 4 ) 2 (2.17) 

consisting of three linear equations and only one quadratic equation. Thus we can 
first solve the linear equations for the coordinates W 1 , W 2 , W 3 in terms of W°, which 
then can be determined by the last quadratic equation. As a result, the coordinates 
of the event P are fully determined in terms of the satellite proper times r and the 
known parameters characterizing their world lines. 

Consider an example in which one satellite is at rest at the origin O and the other 
three move along the three spatial axes. Then the 4-velocities are 

Ui = coshai<9 t + smhaid x , 
JJ% = cosha2<9t + sinha!2<9y , 
Us = cosha3(9 t + sinha3(9 z , 

Ui = dt, (2.18) 
where a,, i = 1, 2, 3, are the rapidities. The system (|2.17|) then reduces to 
= A*!^ — W i + <^ 1 1 i = 1, 2, 3 

= -(W°- t 4 ) 2 + S ij W i W i , (2.19) 
where the notation 

r 4 (t 4 ) 2 -(tM 2 

A^cotha,-— — , & = K . ) ' (2.20) 

t % sinh on 2t 1 sinh on 

has been introduced. The solution of Eq. (|2. 19[) is straightforward 

W l =A i W a + & , (2.21) 

while W° satisfies the quadratic equation 

a(W ) 2 + bW° + c = (2.22) 

with coefficients 

a = 1 - %A J A J , b = -2 (r 4 + % A l &) , c = (r 4 ) 2 - , (2.23) 
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which we will assume to be all nonzero hereafter, plus the additional conditions 
$A < ensuring that P is in the strict causal future of Pa, as stated above. 
The components of the frame vectors (12.12)) turn out to be given by 



d 



dr A L s 2aW° + b 



d l c 



dr A \a 



( 9 



\dr A 



% [Q° A + A>U] > Qa = + — r , (2.24) 



w° 


d 




m 


dr A 






W° 4 


d<S> 1 


dr A 


dr A 



so that the components gAB of the transformed metric follow easily 

9ab = -aUSs + Sij [QaQb + Qa^Zb + ^UQb\ ■ ( 2 - 25 ) 
From Eqs. (|2.24p we see that the quantities £4 are fractional linear functions of W° 
whereas the Q\ are simply linear functions of W°. As a consequence, each of the 
metric coefficients can be cast in the form of a fractional linear function of W°, i.e., 

ciabW + b AB 

$ AB = wo 1 j — ' 2 ' 26 ) 

cabW { > + d AB 

where the coefficients ciab, oab, cab and dAB are four independent functions of the 
emission coordinates t a and of the kincmatical parameters of the satellites. Note that 
it is easy to show that this is true in general, not only for our particular choice (|2.18D 
of satellite motion. 

The same approach we have outlined above will be applied in the next section 
to the more physically interesting case of a metric describing the homogeneous 
gravitational field of the Earth. 

3. Fermi vs emission coordinates 

Consider a generic spacetime metric and introduce a Fermi coordinate system 
(T, X, Y, Z) in some neighborhood of an accelerated world line with (constant) 
acceleration A; the spatial coordinates X, Y, Z are associated with three Fermi- Walker 
dragged axes along the world line while T measures proper time along the world line 
at the origin of the spatial coordinates. Up to terms linear in the spatial coordinates, 
one has (see Eq. (6.18) of Ref. [T3]) 

ds 2 = (r ]Q p + 2AXS J p )dX a dX (3 = -(1 - 2AX)dT 2 + dX 2 +dY 2 + dZ 2 +0(2) ,(3.1) 

valid within a world tube region of radius 1/ A so that \AX\ <C 1 is the condition for 
this approximation to be correct. 

Let the gravitational field of the Earth be represented by the exterior 
Schwarzschild solution, whose metric written in standard coordinates (t, r, 9, 4>) is given 

by 



2M\ , , / 2M\ 1 , , 9/ ~ 
1 \dt 2 + 1 dr 2 + r 2 d9 2 



ds 2 = - ^1 - —J dt 2 + I 1 - —J dr 2 + r 2 {d6 2 + sin 2 6»d</> 2 ) . (3.2) 

The case of Schwarzschild spacetime cannot be treated explicitly, since the geodesies 
are not known in closed analytic form, so that an exact analytic expression for the 
world function cannot be obtained. 

We are interested in casting the metric (|3.2|) in the form ()3.1[) about the world 
line of an accelerated observer at rest at a fixed position on the equatorial plane. The 
map between Schwarzschild coordinates (t, r, 9, (f>) and Fermi coordinates (T, X, Y, Z) 
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Figure 1. The satellite configuration in flat spacetime is schematically shown 
in Fig. (a) by suppressing two spatial directions. The past light cone at a generic 
spacetime point P cuts the two satellite world lines Si, S2 at points Pi and P2, 
respectively. Fig. (b) represents instead the construction of emission coordinates 
in a 1 + 1-dimensional spacetime diagram. 



associated with an observer at rest at r = ro, = 8q = tt/2, 4> — 0o has been first 
derived by Leaute and Linet [2] up to second order corrections in the spatial Fermi 
coordinates llj 
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T 
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1 - 
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(Y 2 + Z 2 ) 
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XY 



1/2 



xz 



0(3) 
-0(3) 



ro 

1 / 2M S 
r o V r Q J 

with the uniform acceleration A entering Eq. (|3.ip given by 

2M^ " 1/2 



^0 



(3.3) 



(3.4) 



I This result has been generalized later to the case of a static observer located at any point on the 
equatorial plane of the Kerr spacetime and to any uniformly rotating circular equatorial orbit by 
Bini, Geralico and Jantzen 1151. 
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The above relations allow one to easily find out the explicit tranformation between 
Schwarzschild coordinates and emission coordinates, once the mapping expressing 
general Fermi coordinates and emission coordinates is constructed. Our results thus 
apply in the case of the Schwarzschild and Kerr spacetimes as well. Furthermore, the 
transformation (|3.3[) can be easily improved by including higher order terms in the 
approximation [15] . making the description even more accurate. 

Let us consider now a set of four satellites whose world lines are chosen in such a 
way that they reduce to the flat spacetime configuration (|2.18| in the limit of vanishing 
acceleration parameter A. All relevant quantities are evaluated in the Appendix. The 
four velocities are given by (see Eq. (|A.12[) ) 

U\ = coshai^T + sinhai<9x + At 1 coshai (2sinhai<9T + coshai<9x) , 
U2 = cosh + sinh a^cV + At 2 cosh 2 cxidx , 
U3 = cosha^dr + sinha3<9z + At 3 cosh 2 a 3 dx , 

Ua=D t + At 4 8 x (3.5) 

to first order in A. By integrating these equations with respect to each proper time 
one easily gets the corresponding world lines 



x 1 


= Si- 


- Air 1 ) 2 coshai ^sinhai^T + ^ coshai<9x 


x 2 


= s 2 - 


- ^(r 2 ) 2 cosh 2 a 2 9x , 


x 3 


= s 3 - 


- ^A(t 3 ) 2 cosh 2 a 3 d x , 


A 4 




- lA^fdx , 



(3.6) 

where the zeroth order quantities are given by Eqs. (|2.4[) and (|2 . 1 8|) . 

Consider then a generic spacetime point P with coordinates X a and a photon 
emitted at the generic point Pa with coordinates X°^ along the world line of the AX\\ 
satellite. The equations of null geodesies are listed in the Appendix (see Eq. (|A.8[) ~). 
Let A be the value of the affine parameter which corresponds to the meeting point P 
according to 

X a ^JC a X + X'X, (3.7) 
where the null vector /C is given by Eq. (|A.10|) , implying that 

X% - X a = -K a l . (3.8) 

Taking the norm of both sides (with the metric components evaluated at P) leads to 

- (t A ) 2 + r\ a pX a XP - 2T A j la pUZXP + AX^X 1 + t a (U\X° - 2C/^A 1 )] = (3.9) 

to first order in A, which generalizes the corresponding Eq. (|2.8[) valid in the case of 
flat spacetime. 

• Emission vs Fermi coordinates 
Searching for solutions of the form 

t a =t a ^+At a ^, (3.10) 

where t a ^ is given by Eq. (|2.9|) with W — > X, we get 

_ A(1) A» X°X 1 +r^(")(^A°-2^A 1 ) 

2 TM0) +r]af3U 2XP ■ {6 - LL) 
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The inverse of the tranformed metric at P follows easily 

AB _ a/3_^_£^f_ _ AB(0) , a AB(l) 

9 ~ 9 dX<*dXP~ 9 +Ag 

to first order in A, where g AB ^ is given by Eq. (|2. 10(1 with W — » X and t a 
and 



yAB{X) _ 



a/3 



- 2X 



1 8t a ^ 8t b ^ 



8X° 8X° 



(3.12) 

> T M°) 

(3.13) 



V dX a dXl 3 dX a dX? 

It is easy to show by a direct calculation that the condition g AA = is preserved. 
The vanishing of the diagonal components of the contravariant spacetime metric once 
written in emission coordinates is actually a general property of real null dual frames 
[TT1 [T2"] , which can be easily explained because the t a = const hypersurfaces are 
tangent to the light cone by construction. 

We are left expressing the spacetime coordinates in terms of the emission 
coordinates. The condition ensuring that the past light cone of P cuts the emitter 
world lines writes as 



il(X A ,X) = 0, 
where the world function is given by 



X° A <X° 



(3.14) 



Q(X A , X) ~ \ [r, aP + A{X\ + X 1 )^] {XI - X a ){X f i - X?) 
= Qn,t(X Al X) + \A{X\ + X 1 )^ - X ) 2 



(3.15) 



to first order in the acceleration parameter A (see Eq. (IA.11[0 . 
• Fermi vs emission coordinates 

Eq. (|3.14|) gives rise to a set of four equations for the coordinates of P. We look for 
solutions of such a system to first order in A, i.e., 

X^W + Aw, (3.16) 

where the solution for W is given by Eqs. (12.2ip - (|2.23j) . The zeroth order equations 
(I3.14p are obviously identically satisfied by the flat spacetime solution W. The 
remaining set of equations for the first order quantities w a is given by 

1. 



= w v A 



Al 



w 



W 



W° - W l A 



lAl 



= w"A z - w z - W^W'A 1 , 
= w°A 3 - w 3 - W°W 1 A 3 , 

= - 2w"(W° - r 4 ) + 25 l3 w l W 3 + W°W 1 {W° - 2t 4 ) , (3.17) 

where the quantities A* are given by Eq. (|2.20| . The corresponding solution turns 
out to be 



w 



W v W l 



w 



= w 3 



The components g A B of the transformed metric turn out to be given by 



dX a dX? 



9ab = g a /3- 



9% + 0(2) , 



(3.18) 



(3.19) 



q t a q t b 

where the zeroth order metric g^ B is given by Eq. ()2.25p and the first order metric 
9ab vanishes. According to the terminology of Ref. [11] the components of the metric 
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tensor are thus "complete" observables, since they are completely determined by any 
given set of four emission coordinates. The latter are instead "partial" observables, 
since they are directly measured quantities. This is true if the spacetime metric, 
i.e., the gravity field, is exactly known (in any coordinate system), implying that 
the system of satellites would constitute an ideal positioning system. In practice the 
spacetime metric is not exactly known, and the satellite system itself has to be used to 
infer it, as discussed by Coll and collaborators [16]. The constellation of satellites can 
thus serve for both positioning and measuring the spacetime metric by equipping the 
satellites with an accelerometer (measuring deviations from geodesic motion) and a 
gradiometer (measuring the strength of the gravitational field). By taking advantage 
from this additional information on the metric an optimization procedure has been 
developed in [16] to obtain the "best observational gravitational field" acting on the 
constellation. 

4. Concluding remarks 

We have considered the metric associated with a generic uniformly accelerated observer 
in any spacetime close enough to the world line of the observer himself. In particular, 
this can be taken as the metric describing the homogeneous gravitational field of the 
Earth. We have expressed this metric in terms of the so called emission coordinates, 
i.e., the four proper times measured along the (timelike) geodesic world lines of four 
satellites, generalizing previous results valid for flat spacetime. The present analysis 
has been mostly motivated by the relevance of using emission coordinates in the 
definition of a relativistic positioning system around the Earth. 

We have considered a particular (symmetric) configuration of satellite motion 
allowing certain simplifications of otherwise more involved formulas. However, our 
results can be easily generalized to arbitrary configurations of satellites corresponding 
to more realistic situations. In fact, the resulting metric and all the possible usages 
associated with it can only be implemented numerically in any case. In this respect, 
we have thus provided an algorithm to construct an emission coordinate system at 
the disposal of a user in the close vicinity of the Earth's surface which also takes into 
account the acceleration due to the Earth's gravity. 
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Appendix A. Geodesies of the Fermi metric 

We give here the general form of both timelike and null geodesies of the Fermi metric 
(|3.ip as well as the expression of the world function, the latter of which is not given 
in the literature. 

The timelike geodesies can be written in the form 



X a ~ S a + As a , 
to first order in the acceleration parameter A, or explicitly 



(A.l) 



T(r) = Cr + T + ACt{P x t + X ) = S° + As 
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X(r) 

Y(r) 
Z(t) 



P x r 



P T - 

P z r- 



1 



- X + -AC 2 t 2 = S X + As 1 
Yq = s 2 , 

Zn = S 3 , 



where 



C= [\ + {P x f + {P Y ) 2 + {p z ) 2 } 112 

The zeroth order quantities correspond to those of Eq. 
S a = U a r - 



where now 



Furthermore 



5 a 
' 



U = Cd T + P x d x + P Y d Y + P Z d z 



= Ci 



(P x t + X )5S + -Ct5? 



(A.2) 
(A.3) 
(A.4) 
(A.5) 

(A.6) 



so that the unit vector tangent to the timelike geodesic world lines turns out to be 



U~U + Au = U + AC[{2P x t 4 

Similarly, the null geodesies are given by 

T(A) = EX + T(0) + AEX[K X X 
1 



X )d T + Crd x ] 



x< 



(A.7) 



X(X) = K X X + X(0) + -AE 2 X 2 , 

Y(X) =K Y X + Y(0) , 
Z(X) = K Z X + Z(0) , 
where A is an affine parameter and 

E = ± [(K x ) 2 + (K Y ) 2 + (K z ) 2 ] 1/2 
The tangent vector to the photon path is thus given by 

K si K + Ak = Ed T + K x dx + K Y dY + K z dz + AE[(2K x X + Xo)d T + EXdx} .(A.10) 

Finally, with these explicit expressions of the geodesies and using the definition 
(|1.2|) , it is easy to obtain the form of the world function 

1 



(A.8) 
(A.9) 



n(x A , X B ) ~ - [ VaP + A(X\ + X B )6°J°] (X2 - X%){X 



X 



q^Xa, x B ) + \a(x\ + x B )(x° A - x 



\2 
B> 



(All) 



to first order in the acceleration parameter A, where Xa and Xb are two generic 
spacetime points connected by a geodesic path. 

Eq. (|A.7[) is quite general. For our purposes we need to specify a set of four 
satellite world lines that reduce to the flat spacetime configuration (|2.18[) in the limit of 
vanishing acceleration parameter A. Their four velocities correspond to the geodesies 
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(|A.2[) , all starting from the origin at each value of their proper times set equal to zero 
(i.e., %=0): 

U x = ^l + (P x ) 2 d T + P x d x + Ay/l + {px^r 1 (^2P x d T + y/l + (P x ) 2 d x ^j , 

U 2 = yjl + (P Y ) 2 d T + P Y dy + A[\ + {P Y ) 2 ]T 2 d X , 

= ^l + (P z ) 2 d T + P Z d z + A[l + (P z ) 2 }t 3 8x , 
U 4 = d T + Ar A d x , (A. 12) 

which give Eq. (|3.5[) after introducing the rapidity parametrization 

P x = sinhai , P Y = sinha 2 , P z = sinha 3 . (A. 13) 
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